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ABSTRACT 

After about 18 years of steadily spinning down, the accretion-powered pulsar 4U 1626-67, 
experienced a new torque reversal at the beginning of 2008. For the present study we have used 
all available Fermi/GBM data since its launch in 2008 June 11 and over 5 yr of hard X-ray 
Swift/B AT observations (starting from 2004 October up to the present time). From 2004 up to 
the end of 2007 the spin-down rate averaged at a mean rate of j> = —4.8 x 10^^"^ Hzs^^ until 
the torque reversal reported here. This second detected torque reversal was centered near MJD 
54500 (2008 Feb 4) and it lasted approximately 150 days. During the reversal the source also 
underwent an increase in flux by a fraction of ^2.5. Since then it has been following a steady 
spin-up at a mean rate of ~ = 4 x 10~^^Hzs^^. We present a detailed long-term timing 
analysis of this source and a long term spectral hardness ratio study in order to see whether 
there are spectral changes around this new observed torque reversal. 

Subject headings: accretion, accretion disks — binaries: close — pulsars: individual (4U 1626-67 ) — 
stars: neutron — X-rays: stars 



1. INTRODUCTION 

The accreting-pow ered pulsar 4U1626- 67 was 
discovered by Uhuru ( Giacconi et al.l 1972 ) . This 
low mass X-ray binary (LMXB) consists of a 
7.66 s X-ray pulsar accreting from an extremely 
low mass co mpan ion (0.04 M© for i = 18°) 



secondary revealed (jMiddleditch et al.l Il98ll) the 
42 mi n orbital period, confirmed bv IChakrabartv 
(|l998l) . Most likely the jinary system contains a 



( Levine et al. 19881 ). Although orbital motion has 
never been detected in the X-ray data, pulsed 
optical emission reprocessed on the surface of the 



hydrogen-deple ted secondary to reach such a shor t 
orbital period ( Paczvnski and Sienkiewic3 198lf ). 
The faint optical counterpart (KZ TrA, V~17.5) 
has a strong UV excess and high opt ical pulse frac- 
tion (jMcClintock et al.lll977l ll980f) . A persistent 
48 mHz quasi-periodic oscillation ( QPO) has been 



detected in the X- 



ray 
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emission (IShmoda et al 



Il990l: iKom mers et al.l[l998l ). In 2008 iKaur et al 



(,20081 claimed that the QPO frequency evolu- 
tion during the previous 22 years changed from 
a positive to a negative trend, somewhat coinci- 
dent wi th the June 1990 torq ue reversal in this 
source. Orlandini et al. ( 1998[ ) inferred a neutron 
star magnetic field in the range (2.4-6.3) x 10^^ G, 
based on the idea that the quasiperiodic oscillation 
frequency is due to the beating between the pulse 
frequency and the Keplerian motion at the magne- 
tospheric radius. To compute this magnetic field 
range a source distance of 5-13 kpc was assumed. 
This distance range was obtained from measure- 
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ments of the optical and X-ray fluxes, assuming 
that LMXB ac cretion disks have a very effective 
X-ray albedo (jChakrabartvl ll998L and references 
therein). Using a ~37 keV absorption cyclotron 
featur e found in the 0.1-200 keV BeppoSAX spec- 
trum. lOrlandini et al. ( 19981 ) obtained a value for 
the magnetic field of 3.2(1 + z)x 10^^ G, where z is 
the gravitational redshift. The X-ray broad-band 
continuum was fitted with a low-energy absorp- 
tion, a blackbody, a power law and a high energy 
cutoff. 
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Fig. 1. — Pulse frequency history of 4U 1626-67, 
showing all available historical data from 1997 to 
2003. The 1990 June (-MJD 48000) torque rever- 
sal is clearly seen.ff 

For mor e than a decade after t he discovery of 
pulsations (jRappaport et al.lll977l ) the source un- 
derwent steady spi n-up at a mean rate of ~ v = 
8.5 X 10-13 Hz s-i (jChakrabartv et al.lll997h (see 
Fig[T]). Monitoring of the source by the Burst and 
Transient Source Experiment (BATSE) on board 
the Compton Gamma Ray Observatory (CGRO) 
starting in April 1991, found the pulsar spinning 
down, i mplying a changed sign in the accretio n 



torque (|Wilson et al.Ml993l : iBildsten et al.l 119941 ). 



During the 7 years after the first torque rever- 
sal, the pulsar spun -down at a rate of 
-7.2 X 10-13 



H s-i (|Chakrabartv et al.lll997n 



Torque reversals have also been observed in 
the accreting pulsar systems Her X-1, Cen X-3, 
GX 1+4, OAO 1657-415, Vela X-1 (IBildsten et al 



19971) and 4U1907-h09 ([in^eEaD Eooi). The 

systems most similar to 4U 1626-67 are Her X-1 
and Cen X-3, where a persistent accretion disk 
is know to be present. Her X-1 like 4U 1626- 
67 has a Roche-lobe filling low mass companion, 
where matter from the companio n flows from the 



LI p oint into an accretion disk (jRevnolds et al 



19971 ). Cen X-3 has a high mass companion which 



almost overflows its Roche-lobe with a focused 
wind from companion flowin g into the accretion 
disk (jPav and Stevens! [1993I) . Studies of accre- 
tion torque in Her X-1 are made difficult by 
the obstruction caused by the warped accretion 
disk which is viewed from near edge-on, result- 
ing in attenuation or comple te eclipsing of the 



pulsar (jPetterson et al.l 119911 ). make fiux mea 



surements unreliable probes of the mass accretion 
rate. In Cen X-3 the dense wind results in a large 
column density, with source occasi onally being 
comp letely obscured. In GX 1-1-4 (iHinkle et a.1 



(jPupree et al 



20061) . OAO 1 657-415 dChakrabartv et all \20o¥. 



4U1907-H09 (iFritz et all l2006h and~\^la X-1 



19801) the optical companion under- 



fills its Roche lobe and accretion proceeds from a 
wind, with perhaps the transient formation of an 
accretion disk if the accreting angular momentum 
is large enough. Vela X-1 is the prototype for a 
super-giant wind fed system, where accretion pro- 
ceeds by direct capture from the wind. In this 
case the transfer of angular momentum is very 
inefficient. Detailed timing analysis show that the 
frequency history of Vela X-1 is consistent with a 
random walk, or equivalently the power-spectrum 
of the torque is white noise (jDeeter et al. 
However, rare short high-flux states have been 
obse rved where transient disk accretion may oc- 
cur ( Krivonos et al. I l2003h . It is likely transien; 



accretion disks occur in GX 1-1-4, OAO 1657-415, 
4U 1907-1-09 and play an essential part there fre- 
quency histories. Understanding the torque be- 
havior of these systems will require characteri- 
zation of the companions wind, and detecting the 
pretense or absence of an accretion disk in addition 
to understanding the flow of angular momentum 
between that disk and the neutron star. 4U 1626- 
67 in contrast provides a cleaner laboratory for 
investigating the flow of angular momentum. 

We present a long term timing analysis using all 
the available Fermi/ GBM data since its launch 
in 2008 June 11 and over 5 yr of hard X-ray 
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Swift/BAT data from 2004 up to 2009. Spectral 
analysis was also carried out for 4U 1626-67, in 
order to see changes during this second detected 
torque reversal that could help us to better un- 
derstand the physical mechanisms involved in this 
process. 

2. Fermt/GBM 

2.1. OBSERVATIONS 

Since 2008 June 11 4U 1626-67 has been con- 
tinuously monit ored by the G amma-ray Burst 
Monitor fGBM) (|Meeganl I2009L submitted), on 
board the Fermi observatory. Timing analysis 
was carried out with GBM CTIME data, with 8 
channel spectra every 0.256 seconds. The total 
exposure time was ~13.75 Ms. 

The GBM is an all-sky instrument sensitive to 
X-rays and gamma rays with energies between ~8 
keV and ~40 MeV. GBM includes 12 Sodium Io- 
dide (Nal) scintillation detectors and 2 Bismuth 
Germanate (BGO) scintillation detectors. The 
Nal detectors cover the lower part of the energy 
range, from 8 keV to about 1 MeV. The BGO de- 
tectors cover the energy range of ~150 keV to ~40 
MeV. Only data from the Nal detectors were used 
in the analysis presented in this paper. 

2.2. TIMING ANALYSIS AND RESULTS 

The analysis of the Fermi/GBM data is com- 
plicated by Fermi's continuously changing orien- 
tation. All intervals of CTIME data from the 
12 Nal detectors are selected for analysis where 
the high voltage is on, excluding those containing 
high voltage transients, phosphorescence events, 
rapid spacecraft slews. South Atlantic Anomaly 
induced transients, electron precipitation events 
and gamma-ray bursts. Source pulses are then 
separated from the background by fitting the 
rates in all detectors with a background model, 
and subtracting the best fit model. This model 
includes bright sources and their changing detec- 
tor responses, including Earth occultation steps, 
and quadratic spline functions which account 
for the remaining long-term background trends. 
The spline models have statistical constraints on 
the changes in second derivative between spline 
segments to control the model stiffness. These 
fits are made jointly across detectors (with com- 
mon bright source fluxes) but separately for each 
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Fig. 2. — Top panel. Fermi/GBM pulse fre- 
quency measurements of 4U 1626-67 since 2008 
August. A change in the sign of the torque was 
found after 18 years of the source spinning down. 
Bottom panel. Swift/BAT pulse frequency history 
covering this second reversal (from 2004 Oct to 
the present time). Error bars are smaller than the 
plotted symbols. 



channel of the CTIME data. Then we combine 
the residuals over detectors with time dependent 
weights which are proportional to the predicted 
(phase averaged) count rates from the pulsar. 
Short intervals (~300s) of these combined resid- 
uals are then fit with a constant plus a Fourier 
expansion to determine a pulse profile. The pro- 
files are divided into six day intervals and the pulse 
frequency and mean profile determined in each in- 
terval with a search of pulse frequ ency for the 
maximum of the Y„ (n=2) statistic (jFinger et al 
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Fig. 3. — Top panel. Swift/BAT spin-up rate his- 
tory of 4U 1626-67 . Middle panel. Average 15- 
50 keV BAT count rate vs. time. Error bars are 
smaller than the plotted symbols. Bottom panel. 
BAT count rate vs. spin-up rate for all the pe- 
riod (circles). A correlation pattern is observed 
specially during the torque reversal (only square 
symbols). 



19991) . Y„ was formulated for finding a pulse 
frequency from a series of pulse profiles, each 
represented by a finite Fourier expansion, and 
accounts for possible frequency dependent non- 
Poisson noise. 

Our monitoring of 4U 1626-67 with Fermi/GBM 
starting in 2008 August, found the pulsar again 
spinning-up rather than spinning-down. Fig [2] 
shows the pulse frequency history using data from 
this monitoring. 4U 1626-67 seems to be increas- 



ing in !>. Follow-up Fermi/GBM observations con- 
firm that the pulsar it is currently spinning-up at 
a mean rate of ~ !>=4xl0'^^^ Hz s^^. 

3. Swift/BAT 

3.1. OBSERVATIONS 



Th e Swift Gamma-ray mission (jGehrels et al 
20041 ) was launched on 2004 November 20. The 



hard X-ray (15-150 keV) Burst Alert Telescope 
(BAT) on board Swift monitors the entire sky 
searching mainly for GRBs. While searching for 
bursts, BAT also accumulates a hard X-ray sur- 
vey of the entire sky covering 2 sr at any par- 
ticular time. Therefore it produces continuous 
streams of rate data. The rate data include quad- 
rant rates (1.6 sec sampling; four energy bands; 
four separate spatial quadrants), not background 
subtracted. For the present study we have ana- 
lyzed more than 4 years of BAT quadrant rates ob- 
servations when 4U 1626-67 was visible (total ex- 
posure time ~ 13 Ms). For the hardness ratio anal- 
ysis we used count rates from the Swift/BAT tran- 
sient monitor results provided by the Swift/BAT 
teairl^. 

3.2. TIMING ANALYSIS AND RE- 
SULTS 

A similar procedure was followed for the 
Swift/BAT quadrant rates timing analysis. Ini- 
tial good time interval (GTI) files are obtained 
using the maketime ftool (hcasoft-6.6.1 jl. Then a 
filtered version of the quadrant rates is obtained, 
rejecting those times when the source is below 
the horizon, to then finally be barycentered us- 
ing the ftool barycorr. Data are then inspected 
and cleaned as in the previous section. With the 
ftool hatmasktgimg the pixel exposure fraction for 
each quadrant is computed for the center of each 
(refined) GTI interval. Pulse profiles for each 
good GTI interval are computed. First the rates 
for each quadrant are fit to a quadratic -|-Fourier 
expansion. Then the Fourier coefficients are com- 
bined using the quadrant exposures to produce 
mean profiles (with units of counts s""'^ cm"'^). In 
a final stage, the Y„ (n=2) statistic is again used 
in intervals of 35 days and a frequency search for 



^ http: / /heasarc. gsfc.nasa.gov/docs / swift / results /transients 
^http: / /heasarc. gsfc.nasa.gov/docs/software. html 
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Table 1: RXTE/PCA and HEXTE CONTINUUM SPECTRAL FITS^ 
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pulsations is carried out. The spin rates were com- 
puted by fitting a linear function to the frequen- 
cies, which were divided into 21 time intervals. 

Swift/BAT observations allowed us to cover the 
evolution of this second torque reversal. We found 
that the pulsar spun-down at a mean rate of 
^ i/ = —4.8 X 10~^'^Hzs~^ until the source re- 
versed torque. Fig [2] shows that the transition 
took place at around MJD 54500 (2008 Feb 04) 
and lasted approximately 150 days. In the bot- 
tom panel of Fig [3] we can see that there is a 
strong correlation between the Swift/BAT count 
rate and the spin-up rate especially during the 
reversal. It is not known if this occurred in the 
1990 torque reversal because of the scarce number 
of observations. In order to see any other possi- 
ble change, we created pulse profiles in the 15-50 
keV band (two harmonics were selected) . We have 
not observed any significant change in pulse shape, 
not even during the reversal. They all are single- 
peaked and sometimes n ot entirely symmetr ic as 
was recently reported by iKrauss et al. ( 2007 ) . 



4. RXTE 

4.1. OBSERVATIONS 

The Ross i X-ra y Timing Explorer {RXTE) 



(|Bradt et al.l Il993f ) carries 3 instruments on 



board. The Propor tional Counter Array (PCA) 
( Jahoda et al. 19961 ) is sensitive from 2-60 keV. 
The High Energy X-ray Timing Experiment 



(HEXTE) (jOruber et all Il996l ) extends the X- 
ray sensitivity up to 200 keV. Monitoring the 
long-term behavior of some of the brig htest X-ray 



sourc es, the All Sky Monitor (ASM) (Levine et al 



19961 ) scans most of the sky every 1.5 hours at 2- 
10 keV. Two RXTE/PCA observations from 2008 
March 5 and 13 were used (ID 93431-01-01-00 
and 93431-01-02-00; 7.174 ksec). For spectral 
analysis we selected PCA Standard-2 data which 
contains 129-channel range spectra taken every 16 
seconds and HEXTE Standard Modes (Archive) 
data which contains Spectral Bin (64-bin spectra 
produced every 16s). For the long-term hard- 
ness ratio analysis we used the ASM daily flux 
averages in the 1.5-12 keV energy range from the 
HEASARC archivel. 

4.2. SPECTRAL ANALYSIS AND RE- 
SULTS 

RXTE/PCA (2.5-20 keV) and HEXTE (18- 
100 keV) spectra were fitt ed in XSPEC 11.3 2 
with two models used by IPravdo et al. (Il979l ). 
Using these models allows us to com pare our spec- 



tral s t udy with previous wo r ks bv (iPravdo et al 



1979t lOrlandini et all Il998l; iKrauss et al.l 120071: 



Jain et a l. 2009) and update the long-term X- 
ray flux history of 4U 1626-67 relative to the flux 
measured by HEAO 1 (Chakrabarty et al. (1997); 
Krauss et al. (2007)) . The first model includes 



^http: / /hcasarc. gsfc.nasa.gov/docs/archive. html 
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a low-energy absorption, a blackbody component, 
a power law and a high-energy cutoff at ~20 keV 
(WABS (GAUSS-l-BBODY-fPOWLAW) HIGHE- 
CUT). A broad line near 6.5 keV significantly im- 
proves the present fit and indicates the presence 
of an iron line, also suggested by IPravdo et alj 
(|l979l ) in their ( 0.7-100 keV) spectral analysis 
of this source. The column density of cool ma- 
terial in the line of sight was fixed in our study 
since it could not be constrained. A value of 
1.3 X 10^^ cm~^ was selected from Krauss et al.l 
(|2007n . The spectral parameters obtained are 
shown in Table 1 . Reprocessing of photons either 
at the base of the accretion column or the inner 
edge of the accretion disk might explain the 
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Fig. 4. — Top panel. Long term hardness ratio 
analysis of 4U 1626-67. Swift/BAT count rates 
(15-50 keV) were selected as the hard band and 
RXTE/ ASM count rate (1.5-12 keV) as the soft 
band. Bottom panel. Hardness-intensity dia- 
gram. During the reversal a transition from hard 
to soft is seen. 



origin on the blackbody component. From the 
blackbody flux and temperature we obtain an 
emission area of 9 x 10"'^^ cm~^ (assuming a source 
distance of 10 kpc), which is inconsistent with ei- 
ther of those possibilities. We flt in addition the 
same model with a bremsstrahlung instead of a 
blackbody component , obtaining a compatible flt. 
Table [T] summarizes the spectral parameters ob- 
tained. Fluxes for the first model in the 2-10 keV, 
2-20 keV and 2-50 keV bands are 3.24(2) x 10"!°, 
7.003(11) X 10-1° and 9.98(2) x 10"!° erg cm-^s"!, 
for the first RXTE observation (similar values for 
the second one). 

5. HARDNESS RATIO ANALYSIS 

Fig m (top panel) shows the hardness ratio 
(HR) analysis carried out for this source. The 
HR was defined as the ratio 15-50keV/1.5-12keV 
(BAT/ ASM). To reduce large uncertainties the 
light curves were rebinned and then the HR were 
computed. Like all hardness ratios, these are in- 
strumentally dependent. We can see that there is 
a smooth hardening evolution of the source before 
the reversal, although from this figure we cannot 
observe any dramatic change a posteriori. 

In the bottom panel of Fig |4] we have replaced 
the intensity with the BAT count rate in order to 
perform Hardness-intensity diagram (HID). This 
allow us to study the long-term spectral variability 
of 4U 1626-67 , including the transition, since the 
2 RXTE observations do not provide us any direct 
comparison between before and after the torque 
reversal. From that figure we can see that there 
is a transition from hard to soft during this new 
reversal of 4U 1626-67. 

6. DISCUSSION 

The discovery of the torque reversal in 2008 (see 
FigO has shown the spin behavior of 4U 1626-67 
to be more complicated than previously thought. 
While the continuous decrease of total X-ray flux 
during almost two decades was expected to br ing 



the source into quiescence (jKrauss et al.l 120071 ) . a 



new spin reversal with a rapid increase of the flux 
occurred. This inconsistency reopens discussion 
on physical processes governing the spin evolution 
of the pulsar, and raises a question about the na- 
ture of the new torque reversal. 

All previous studies of 4U 1626-67 were focused 
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on modeling the spin-up torque applied to the neu- 
tron star from the accreted material. It was widely 
believed that the spin behavior of the pulsar de- 
pended mainly on variations of the mass accretion 
rate onto the stellar surface and therefore, the rate 
of mass transfer between the system components. 
An attempt to explain the torque reversal in 1990 
(the spin-up to spin-down transition) in terms 
of va riable equilibrium p eriod (Ghosh and Larnb 
I1979I) has been made by IVaughan and Kitamoto 
(|l997l ). The neutron star phase transition in their 
model is associated with a substantial decrease of 
mass accretion rate and, possibly, change of the 



parameters of 4U 1626-67 is (|Lamb et al.lll973l ) 
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Fig. 5. — Pulse frequency history of 4U 1626-67 
from 1997 up to 2009. The 1990 June and the 
2008 Feb reversals are clearly seen. 

structure of the accretion flow beyond the mag- 
netosphere. This scenario is based on some evi- 
dence for a connection between torque and X-ray 
luminosity in 4U 1626-67 and a hardening of the 
energy spectrum after the torque reversal. They 
have emphasized, however, that this approach en- 
counters major difficulties explaining the observed 
behavior of QPOs, which have been interpreted in 
terms o f the magn ctosphcric beat-frequency oscil- 
lations (jLamb et al . 1985). Their analysis suggests 
that either this interpretation of QPOs is incorrect 
or the torque reversal in 1990 occurred as the mag- 
netospheric radius of the neutron star, r^, had in- 
creased beyond its corotation radius, which for the 



GM„ 



\47r2i 



1/3 



6.7 X 10® m^/^Pj{^ cm. (1) 



Here m and P7.7 are the mass and spin period of 
the neutron star in units of 1.4 M© and 7.7 s. A 
dramatic change in the power spectra between the 
last observation of 4U 1626-67 duri ng the spin- 
down phase (2003) (jKaur et al.l 120081 ). and obser- 
vations made soon after the n ew torque r eversa l 
has been recently reported by IJain et al. ( 20091 ) . 
with the 35-48 mHz QPO no longer being present, 
and wide shoulders on the pulse fundamental ap- 
pearing. They claimed that the observed behavior 
of the source cannot be a simple case of increased 
mass transfer rate, but is also a change in the ac- 
cretion flow parameters. 

Analyzing the evolution of the source en- 
ergy spectrum and possible correlation between 
the torque and X-ray luminosity of the pulsar, 
Yi and Vishniac ( 19991 ) proposed a scenario in 
which the torque reversal in 1990 is associated 
with a state transition of the accretion disk to a ge- 
ometrically thick, hot and, possibly, sub-Keplerian 
phase. Following this idea one could associate 
the 2008 torque reversal with an inverse transi- 
tion of the disk into its previous geometrically 
thin Keplerian phase. However, the reason for 
such a transition is rather unclear since the level 
of X-ray flux measured before and even after the 
2008 reversal is smaller than that measured during 
the reversal in 1990. Furthermore, both reversals 
have occurred at almost the same timescale (about 
150 days), which significantly exceeds the dynami- 
cal timescale in the hot disk in which its transition 
to the ground state is expected. A difficulty to fit 
in the transition timescales o bserved in 4U 1626- 
67 ha, s also been mentioned by Wijers and Pringl3 
(Il999l) . who discussed a possibility to explain the 
torque reversals in terms of the warped disk tran- 
sition into a retrograde regime. 

A correlation between the torque applied to the 
neutron star in 4U 1626-67 and X-ray flux of the 
system in the above mentioned models has been 
adopted as one of the basic assumptions. To test 
the validity of this assumption using data derived 
before 1993 was rather complicated. This is il- 
lustrated in Figure El which shows the 4U 1626- 
67 X-ray flux history. We can see all previous 
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flux measurements and two RXTE/PCA recent 
values from the present work (in the 2-20 keV 
band). These values are relative to the flux mea- 
sured in 1978 by HEAP 1 i n the same energ^i 



surea m iMiii py ntjAU i m tne same energy 
band (lOrlandini et alj llQQSt iKrauss et ahl 120071: 
Chakrabartv et al. 19971 and references therein). 
The cross point before the 2008 reversal has been 
inferred by scaling the PCA fluxes according to 
the observed change (2.5 factor) in the Swift/BAT 
rate, since no spectral changes across the transi- 
tion have been observed according to the present 
work. It should be noted, however, that all these 
relative flux values were computed in different en- 
ergy bands, therefore the general decreasing trend 
we see since 1977 might not be the real picture. 
Points could misrepresent the HEAO 1 0.7-60 keV 
flux by a ^ 20% due to spectral changes. More- 
over, rather than a continuous decline (ignoring 
the 1990 point) 4U 1626-67 might present a flat 
behavior until the first reversal, then a sudden 
drop after that and a decreasing tendency until 
the 2008 reversal, in which this source experienced 
a rapid increase of flux. It therefore appears that 
the spin-down phase is the only suitable part of the 
light curve for testing the correlation between the 
torque and mass accretion rate onto the neutron 
star surface. 

As seen from Figure |6l the X-ray flux during 
the spin-down phase has decreased by a factor of 
2. This indicates that the mass accretion rate onto 
the surface of the neutron star, M, and, corre- 
spondingly, the spin-up torque applied to the star 
(jPringlc and Rccs 1972), 



Af(GM„,r^)i/2, 



(2) 



during this phase have also decreased by at least 
the same value. If the spin-down torque applied 
to the neutron star during this time were constant 
one would expect the pulsar to brake harder at 
its fainter state close to the end of the spin-down 
phase. However, observations show the situation 
to be just the opposite. The spin-down rate of the 
neutron star during th is phase has decreased fr om 
\u\ ~ 7xlO-"Hzs-i (jChakrabartv et al.lll997l) to 
5 X lO^^'^ Hz s^^ (see Figure ??), implying that the 
pulsar was braking harder at its brighter stage just 
after the torque reversal in 1990. According to the 
equation governing spin evolution of an accreting 
neutron star, 



this means that the spin-down torque, Ksd, during 
the spin-down phase has been decreasing simulta- 
neously with the spin-up torque but at a higher 
rate and, therefore, the pulsar spin evolution dur- 
ing this time has been governed mainly by varia- 
tions of Ksd rather than Kgu (here / is the moment 
of inertia of the neutron star) . This conclusion se- 
riously challenges the possibility of modeling the 
spin history of 4U 1626-67 solely in terms of vari- 
ations of M, and suggests that the dramatic in- 
crease of X-ray flux observed in 2008 torque re- 
versal may be a consequence rather than a reason 
for this event. 

With the lack of correlation between the X-ray 
flux and the torque applied to the neutron star, 
modeling of the spin-down torque appears to be 
the main target for theoretical studies of the sys- 
tem. Unfortunately, this part of modeling of the 
magneto-rotational evolution of neutron stars re- 
mains so far a work in progress. The canonical pre- 
scription for the spin-down torque ( Lipunov.,199S} ). 



(4) 



27r/i> = 



(3) 



in our case turns out to be rather ineffective. 
The dipole magnetic moment of the neutron star, 
during the spin-down phase obviously remains 
constant and the corotation radius, changes only 
by 0.16% (Arc < 1.1 x 10^ cm) as it follows 
from the observed changes of the pulsar frequency 
(|Az/| ~ 3.3 X lO"**). Hence, any variation of Ksd 
under these conditions proves to be determined 
by the dimensionless parameter fct < 1, which 
in the case of disk accretion is just t he a pa- 
rame ter at the inner radius of the disk (jLipunov 
19921 . and references therein). However, it ap- 
pears to be rather difficult to combine the assump- 
tion about significant variation of turbulence at 
the inner radius of the disk with the extremely 
low level of noise strength observed in 4U 1626-67 
(jChakrabartv et al.lll997l) . Furthermore, different 
rates of variations of Ksd and Ksn during the spin- 
down phase indicates that either fct is a non-linear 
function of a, or the expression ([4]) is oversimpli- 
fied, and in the case of 4U 1626-67 is ineffective. 
As recently shown bv lPerna et al.l ( 2006f ). the spin- 
down torque proves to be a strongly dependent on 
the relative velocity between the disk and magne- 
tosphere and its value significantly varies even for 
a relatively small variations of the mass-transfer 
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rate if the interaction between that disk and the 
magnetosphere occurs at the corotation radius of 
the star. They proposed a model where simulta- 
neous with accretion from a disk onto the neutron 
star some material from near the disk - magneto- 
sphere boundary is ejected and either escapes from 
the system or is recycled back into the accretion 
disk. This results in a hysteresis-type limit cycle 
where slow changes in the accretion rate from the 
companion into the accretion disk can result in 
the rapid change in torque and luminosity. Their 
model predicts, however, that the luminosity af- 
ter a spin-down to spin-up torque reversal would 
be higher than the luminosity after a spin-up to 
spin-down torque reversal, which is the opposite of 
what occurred for 4U 1626-67 for this new reversal. 
Moreover, for 4U 1626-67 they predicted the full 
spin-down/spin-up cycle would take thousands of 
years, again inconsistent with the present observa- 
tions. It should be noted, however, that the cycle 
could be interrupted by a sudden increase of the 
mass transfer rate in the disk. The brightening of 
the pulsar observed during the torque reversal in 
2008 does not allow us to exclude this possibility. 

The above mentioned problems may indicate 
that the torque applied to the star should be 



treat ed in a different way (see, e.g. lRappaport et al 
20041 ). or the interaction between the disk and 
magnetosphere in 4U 1626-67 occurs in a particu- 
lar region, in which small perturbations of pulsar 
parameters may lead to dramatic changes in the 
torque and mass a c cretion rate . As shown by 
Anzer and Boerner ( 1980l 1983 ). the boundary 
between the disk and magnetosphere is subject to 
Kelvin-Helmholtz instability in a region specified 
with the condition 



<<5c 



(5) 



while apart of this region the interchange-type in- 
stabilities of the boundary are suppressed. Here 
5c is the thickness of region where the velocity 
difference between the disk and magnetosphere 
is smaller than the sound velocity of plasma in 
the disk, i.e. \Vi,ir^ ± Jc) - M<(rc)| < V;. Here 
Vk{r) = {GMns/rY^^ is the Keplerian velocity, 
and Vs is the sound speed at the inner radius of 
the disk. The value of 5c is model dependent, but 
eventually is of t he order of the disk th ickness 
(Anzer and Boerner 1983: 



which is comparable with the amplitude of vari- 
ation of the corotation radius during the spin- 
down phase 1990-2008 inferred from our observa- 
tions. This represents some grounds to assume 
that the spin and intensity evolution of 4U 1626- 
67 can be governed by variation of relative posi- 
tion of Tc and Tin in the region specified by con- 
dition Since Arc — (l/2)i5c one could expect 
that the spin- up/spin-down cycle of the pulsar is 
determined mainly by variations of its corotation 
radius, while the magnetospheric radius remains 
almost constant. The recurrent time of the cycle 
in this case would be ~ 15 — 25 year, which is close 
to the observed value. Further analysis of the cor- 
responding scenario is, however, beyond the scope 
of this paper. Here we would like to note only that 
4U 1626-67 represents an exceptional case as an 
accretion-powered pulsar with an extremely small 
torque noise. This indicates that the accretion pic- 
ture in this system may differ significantly from 
that realized in other pulsars. 
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Fig. 6.— The X-ray flux history of 4U 1626-67 rel- 
ative to the flux measured by HE AO 1, in the same 
energy band, from previous works (Chakrabarty 
et al. (1997): circles; Orlandini et al. (1998): 
triangle; Krauss et al. (2007): stars) and two re- 
cent RXTE/PCA observations (unfilled squares) 
in the 2-20 keV band. The cross point is inferred 
from PCA flux and the fractional change in the 
Swift/BAT rate, since no spectral changes during 
the transition have been observed in this work. 
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Finally, the spectral evolution of 4U 1626- 
67 during the torque reversal differs from that 
expected in models which suggest significant 
changes of the accretion flow s tructure in spin- 
up/spin-down transitions (e.g . Yi and Vishniad 
19991 : IWiiers and Prindd \l9m . As seen from 
Figure [H the spectrum becomes the hardest dur- 
ing the reversal and the value of the hardness ratio 
before and after these events does not differ sig- 
nificantly. This indicates that the recent torque 
reversal can be associated with changes of physi- 
cal conditions in the inner part of the disk or/and 
in the region of its interaction with the magne- 
tosphere rather than a significant change of the 
accretion flow geometry. The errors of the obser- 
vations are, however, too large for a justification 
of particular transition model. A more precise 
spectral measurements of the pulsar during its 
next spin- up/spin-down transition is, therefore, 
strongly desired. Assuming the recurrent time of 
the transition to be about 18 years one can sug- 
gest to pay more attention to the pulsar in 2025- 
2028. Since the typical duration of the transition 
is about 150 days, a regular monitoring of the pul- 
sar, frequently enough to provide a spin-up rate 
measurement every two months, would prevent us 
from missing its next torque reversal. 

7. CONCLUSIONS 

We report on a discovery of a new spin-down 
to spin-up torque reversal in 4U 1626-67. It oc- 
curred after about 18 years of the pulsar's steadily 
spinning down and was centered on 2008 Feb 4. 
The transitions was lasted ^^150 days and accom- 
panied by an increase in the Swift/B AT count rate 
of a 2.5 factor ('^150%). The pulsar spectrum was 
harder during the torque transition than before or 
after. A strong correlation between torque and 
luminosity is inferred only during the transition. 
The spin-up and spin-down rates before and af- 
ter the transition were almost identical (^| \ = 
5 X 10""'^^ Hz s~^). However, the pulsar was brak- 
ing harder at the beginning of the spin-down epoch 
in 1990 than at its end in 2008. Furthermore, the 
spin-down rate during this epoch was decreasing 
simultaneously with the decreasing of the source 
X-ray luminosity. Finally, the spin-down to spin- 
up torque reversal in 2008 has occurred at lower 
luminosity as the spin-up to spin-down torque in 
1990. These properties cannot be explained with 



existing models and appear to be a clue for further 
progress in understanding the mechanism govern- 
ing the torque reversals in the accretion -powered 
pulsars. 
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